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ABSTRACT 

A well known property of the 7-ray sources detected by COS-B in the 1970s, by the Compton 
Gamma-ray Observatory in the 1990s and recently by the Fermi observations is the presence of radio 
counterparts, in particular for those associated to extragalactic objects. This observational evidence 
is the basis of the radio-7-ray connection established for the class of active galactic nuclei known as 
blazars. In particular, the main spectral property of the radio counterparts associated with 7-ray 
blazars is that they show a flat spectrum in the GHz frequency range. Our recent analysis dedicated 
to search blazar-like candidates as potential counterparts for the unidentified 7-ray sources (UGSs) 
allowed us to extend the radio-7-ray connection in the MHz regime. We also showed that below 
1 GHz blazars maintain flat radio spectra. Thus on the basis of these new results, we assembled 
a low-frequency radio catalog of flat spectrum sources built by combining the radio observations of 
the Westerbork Northern Sky Survey (WENSS) and of the Westerbork in the southern hemisphere 
(WISH) catalog with those of the NRAO Very Large Array Sky survey (NVSS). This could be used 
in the future to search for new, unknown blazar-like counterparts of the 7-ray sources. First we 
found NVSS counterparts of WSRT radio sources and then we selected flat spectrum radio sources 
according to a new spectral criterion specifically defined for radio observations performed below 1 
GHz. We also described the main properties of the catalog listing 28358 radio sources and their 
logN-logS distributions. Finally a comparison with with the Green Bank 6-cm radio source catalog 
has been performed to investigate the spectral shape of the low-frequency flat spectrum radio sources 
at higher frequencies. 

Subject headings: galaxies: active - quasars: general - surveys - radiation mechanisms: non-thermal 


1. INTRODUCTION 

Since the epoch of the first 7-ray observa¬ 
tions performed by COS-B in the 1970s (e.g., 
iHermsen et al. 1977f ) and by the Compton Gamma-ra y 
Observatory in the 1990s (e.g., [Hartman et al. 19991) . 
a link between the radio and the 7-ray sky was 
discovered. It has been used to associate the high- 
energy so urces wit h their low-energy counterparts 
(e.g. iMattox et al. 1997h . This radio-to-7-ray re¬ 
lation has also been recently highlighted for the 
extragalacti c source s detected by the Fermi mission 
llAtwood et al. 20091 ). In particular, nearly all the 7-ray 
sources associated in the seco nd Fermi Large Area Tele¬ 
scope (LAT) catalog (2FGL; iNolan et al. 2012h and/or 
in the secon d catalo g of active galactic nuclei (AGNs) 
(jAckermann et al. 2011af) detected by the Fermi -LAT 
have a clear radio counterpart. This is the basis of the 
radio-7- ray conn ectio n specific ally discussed for bla zars 
(e.g., iGhirlanda et al. 20101 iMahonv et al. 201(1 
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lAckermann et al. 201 lbf ) , that constitute the 
rarest class of AGNs (e.g., Urrv fe Padovani 19951 
[Massaro et al. 20091 : IMassaro et al. 201 lal ) and the 
largest known population of 7-rav sources (e.g., 
lAbdo etah 201(lh 

Recently we addressed the problem of searching for 
7-ray blazar candidates as counterparts of the unidenti¬ 
fied 7-ray sources (UGSs) adopting a new approach that 
employs the low-frequency radio observations performed 
by the Westerbork Synthesis Radio Telescope (WSRT). 
While performing this investigation we found that the 
radio-y-r ay connectio n of blazars can be extended below 
~1 GHz (IMassaro et al. 2013af) . 

Our analysis was based on the combination of the 
radio observat ions from Westerbork Northern Sky Sur¬ 
vey (WENSS; iRengelink et al. 19971) at 325 MHz with 
those of the NRAO Very Large Array Sky survey 
(NVSS; lCondon et al. 19981) and of the Very Large Array 
Faint I mages of the Radio Sky at Twenty-Ce ntimeters 
(FIRST: iBecker et al. 19951 iWhite et al. 199711 at about 
1.4 GHz. A similar analysis was also performed using the 
Wester bork in the sou thern hemisphe re (WISH ) surve y 
(|De Breuck et al. 2002fl at 352 MHz (jNori et al. 201411 . 
Both of these studies were based on the observational 
evidence that blaz ars a lso sh ow flat radio spectra belo w 
~1 GHz (see also iKovalev 2009at Kovalev et al. 2009bl : 
iPetrov et al. 20131 for recent analyses). 

The flatness of the blazar radio spectra is a 
well known property expected from radio data in 
the GHz frequen c y range (e.g.. Ilvezfc et al. 20021 
iHealev et al. 20071 : iKimball fe Izevic 2008L for recent 
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analyses). This spectral property was also used in 
the past for the as sociations o f y -ray sources since the 
EGRET era (e.g., iMattox et al. 1997}) . However, de- 
spite a small survey of BL Lac objects at 102 MHz 
(|Artvukh fe Vetukhnovskava 198lll . the low radio fre¬ 
quency spectral behavior of blazars was still an unex¬ 
plored region of the electro magnetic spectr um until our 
recent analyses dMassaro et al. 2013at lNori et al. 20141) . 
Using WSRT data at 325 MHz and at 352 MHz as 
well as those of very low-frequency observations of the 
Very L arge Array Low-Frequency Sky SurvejQ (VLSS; 
iCohen et al. 2007? ) at 74 MHz we showed that blazars 
maintain a flat radio spectrum even below ~100 MHz 
and w e extended the radio-y-ray connection below ~1 
GHz dMassaro et al. 201 311) . 

Thus, motivated by these recent results we assembled 
a catalog of low-frequency flat spectrum radio sources us¬ 
ing the combination of both the WENSS and the WISH 
surveys with the NVSS. The main aim of this inves¬ 
tigation is to provide the counterpart, at longer wave¬ 
lengths, of the Combined Radio All-Sky Targeted Eight- 
GHz Survey (CRATES) used to associate F ermi objects 
with blazar-like sources dHealev et al. 20071 ). 

The paper is organized as follows: in Section [2] we 
briefly present the main properties of the low-frequency 
radio survey performed by WSRT and used to carry out 
our investigation (i.e., the WENSS and the WISH). In 
Section[3]we search for the NVSS counterparts of WSRT 
sources. Then in Section |4] we extract the main low- 
frequency catalog of flat spectrum radio sources (LOR- 
CAT) from the combined WSRT-NVSS surveys and we 
discuss on its main properties. Section [5] is devoted 
to the comparison the Green Ba nk 6-crn (GB6) radio 
source catalog (e.g., [Gregory et al. 1996i l to investigate 
the spectral behavior of LORCAT sources at higher fre¬ 
quencies. Finally, Section [6] is dedicated to the summary 
and the conclusions. 

For our numerical results, we use cgs units unless 
stated otherwise Spectral indices, a , are defined by flux 
density, S„ oc v~ a . The WSRT catalogs used to carry out 
our analysis are available from both the HEASARC00 
and the VIZIEB FIFI databases as well as that of the 

NVS330 

2. WESTERBORK LOW-FREQUENCY RADIO SURVEY 

The Westerbork Northern Sky Survey (WENSS) is 
a low-frequency radio survey that covers the northern 
sky above +30° in declination performed at 325 MHz 
to a limi ting flux den sity of ~18 mJy at 5 sigma level 
(iRengelink et al. 199711 . The version of the WENSS cat¬ 
alog used in our analysis was implemented as a combina¬ 
tion of two separate catalogs obtained from the WENSS 
Websit<£3: the WENSS Polar Catalog that comprises 
18186 sources above +72° in declination and the WENSS 
Main Catalog including 211234 objects in the declination 
range between +28° and +76°. 


We also used the Westerbork In the Southern Hemi¬ 
sphere (WISH) catalog[j3 that is the southern extension 
of the WENSS. WISH is a low-frequency (352 MHz) ra¬ 
dio survey covering most of the sky between -26°and - 
9° at 352 MHz to the same limiting flux density of the 
WENSS. It is worth noticing that the Galactic Plane 
region at galactic latitudes |6| <10° are excluded from 
the WISH observations. Due to the very low elevation 
of the observations, the survey has a much lower res¬ 
olution in declination than in right ascension. A cor¬ 
relation with the NVSS shows that the positional ac¬ 
curacy is less constrained in declination than in right 
asce nsio n, but th ere is no significant systematic error 
(see iDe Breuck et al. 200M for more details). Finally, 
we highlight that the WISH catalog contains multiple 
observations of the same source for many objects as 
well as measurements of individual components of multi- 
component sources. 

3. RADIO SPATIAL ASSOCIATIONS 

We adopted the following statistical approach to find 
the radio NVSS counterparts at 1.4 GHz for the sources 
in the WSRT low radio frequency surveys, namely: the 
WENSS and the WISH. 

For each radio source listed in either the WENSS and 
the WISH surveys, we searched for all the NVSS counter¬ 
parts that lie within elliptical regions that corresponds to 
the positional uncertainty at 95% level of confidence (i.e., 
2a). We took into account the uncertainties on both the 
right ascension, a, and the declination, 5, in the WSRT 
and in the NVSS surveys. 

We found that the total number of correspondences is 
225933 out of 268425 radio sources included in either the 
WSRT surveys. We excluded from our analysis all the 
WSRT source with radio analysis flags (i.e., P and Y 
as reported in the WENSS and WISH catalog, respec¬ 
tively, to indicate that there were problems in the model 
fitting for a source) and variability flag in the WISH ob¬ 
servations, all the double matches and all those sources 
labeled as components of a multi-component source (flag 
“C”) in the WSRT catalogs. In addition, for this version 
of the LORCAT catalog, we also excluded from our sam¬ 
ple 2707 multiple matches since their WSRT radio flux 
densities could be due to the emission of several, unre¬ 
solved, NVSS sources so contaminating our estimates of 
the low frequency spectral index. 

We then built 100 mock realizations of the WSRT cat¬ 
alog by shifting each source position in a random direc¬ 
tion of the sky by a fixed length of 1°. This shift adopted 
to create the mock WSRT catalogs were chosen not too 
distant from the original WSRT position and within the 
NVSS footprint so to obtain fake catalogs with a sky dis¬ 
tribution similar to the original WSRT and to perform 
the cross-match with each fake catalog and the NVSS 
taking into account the local density distribution of the 
WSRT radio sources. The total number of WSRT sources 
in each mock realization is also preserved. 

For each mock realization of the WSRT catalog, we 
counted the number of associations with the NVSS occur- 


7 http://lwa.nrl.navy.mil/VLSS/ 

8 WENSS http://heasarc.gsfc.nasa.gov/W3Browse/all/wenss.html 

9 WISH: http://heasarc.gsfc.nasa.gov/W3Browse/all/wish.html . . _ .. . „ . 

10 WENSS: http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=VIII/62'" ln § at an S ular separations R smaller than 300 . Then 

11 WISH: http://cdsarc.u-strasbg.fr/viz-bin/Cat7VIII/69A we computed the mean number A (R) of these mock asso- 

12 http://heasarc.gsfc.nasa.gov/W3Browse/all/nvss.html ciations, averaged over the 100 fake WSRT catalogs and 

13 http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=%20NVSS 

14 http://www.astron.nl/wow/testcode.php?survey—1 15 http://■www.astron.nl/wow/testcode.php7survey=2 
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Fig. 1.— Upper panel) The values of A A (R) (red circles) and 
AN(R) (black squares) as function of the angular separation R. 
Our choice of Rmax is marked by the vertical dashed line. It occurs 
at the first R value for which A A (R) ~ A N(R). Lower panel) The 
probability of having spurious associations p(R) as function of the 
angular separation R. 

verifying that A (R) has a Poissonian distribution. In- 
creasing the radius by AR =5", we also computed the 
difference A A (R) as: 

A X(R) = A(i? + A i?) - X(R ), (1) 

In Figure H] we show the comparison between A N(R) 
and AA(-R). For radii larger than R ma x =95"tlre 
A A (R) curve superimposes that of A N(R) indicating 
that WSRT-NVSS cross-matches could occur by chance 
at angular separations larger than R max ■ Thus we choose 
Rmax as to the maximum angular separation between the 
WSRT and the NVSS position to consider the 1.4GHz 
radio source a reliable counterpart of a WSRT object. 

In addition we calculated the chance probability of spu¬ 
rious associations p(R) as the ratio between the num¬ 
ber of real associations N(R) and the average of those 
found in the mock realizations of the WSRT catalog 
A (R), corre sponding to a value of ~1 0% fo r R = R m „x 
(see e.g.. iMaselli et al~ 201CH IMassaro et al. 20115 

iD’Abrusco et al. 20131 IMassaro et al. 2013d for a sim¬ 
ilar procedure to estimate the probability of spurious as¬ 
sociations). 

We then computed the uncertainties on the WSRT po¬ 
sitions according to the procedure described in Rengelink 
et al . (1997) and combined them with the NVSS ones 
(iCondon et al. 19981) using the following relation: 

&RA,Dec = \J' a RA,Dec(yVSRT) + <J 2 RA Dec (NVSS), 

( 2 ) 

We also defined the angular separation normalized to 
the values of the positional uncertainties m as: 

where Rra and Rdec are the angular separations in 
right ascension and in declination, respectively. In Fig¬ 
ure [2] we show the distributions of the ratio between the 


o> 

-o 


E 


WSRT - NVSS angular separations 




0 0.5 1 1.5 2 2.5 3 

R /a 


Fig. 2.— Left panel) The distributions of ratio between angular 
separations R and the positional uncertainties for the right ascen¬ 
sion (black straight line) and for the declination (red dashed line), 
respectively for the selected 224438 WSRT - NVSS radio sources. 
Right panel) The cumulative distribution of the ratio R/a. The 
dotted blue line marks the 90% limit (see also Section 

angular separation R and the combined positional un¬ 
certainty a for RA and Dec, respectively. In order to 
build the final sample of WSRT - NVSS correspondences 
that will be used to extract the low-frequency radio cat¬ 
alog of flat spectrum sources, we selected only sources 
with m < 3. This WSRT - NVSS final sample lists 
224438 radio sources out of 225933 previously selected. 
We note that we found a potential NVSS counterpart for 
about 85% of the WSRT sources, and since we adopted 
a threshold on m =3, this corresponds to a completeness 
C of about 80%, evaluated according to the relations de¬ 
scribed in Condon et al. (1975). Moreover, this is also 
in agreement with the reliability of our associations, es¬ 
timated via Monte Carlo simulations, occurring at R m ax 
that is of the order to 10% (see Figure [U). 

4. LOW-FREQUENCY RADIO CATALOG OF FLAT 
SPECTRUM SOURCES 

4.1. Radio spectral index distribution at low-frequencies 

For the WSRT-NVSS associations we defined a low- 
frequency radio spectral index: cq ou ,, using the inte¬ 
grated flux densities at 325 MHz from the WENSS and 
those at 352 MHz reported in the WISH, S 325 and S 352 , 
respectively, in combination with the NVSS S 4400 at 1.4 
GHz as: 

c%iow — k± * log(S 1400 /Si ow ) 5 ( 4 ) 

where the ki factor is equal to 1.58 and 1.67 
(i.e., [log{ 1400/325)] -1 and [Zo<7(1400/352)] _1 ) for the 
WENSS and the WISH surveys, respectively, Si ow is 
the flux density at 325 MHz (WENSS) or at 352 MHz 
(WISH) with all flux densities in units of mJy. The un¬ 
certainties on ai ow where computed according to the fol¬ 
lowing relation: 

&low = k 2 ■ \J (cri40o/'S'l40o) 2 + (cq ow /Si OW ) 2 (5) 

where the fc 2 factor is equal to 0.68 and 0.72 (i.e., 

| [/n(1400/325)]- 1 | and | [/n(1400/352)]- 1 |) for the 
WENSS and the WISH surveys, respectively while 0-1400 
and ai ow are the uncertainties on the WSRT and NVSS 
flux densities. 
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Fig. 3. — The distributions of oq ow for three flux limited sub¬ 
samples extracted from the WSRT-NVSS correspondences. Flux 
density cuts are reported in each panel. It is clear how the distri¬ 
bution has a tail for low values of the low-frequency spectral index 
but it does not appear to have a bimodal sh ape as in previous 
radio analyses at higher frequencies (see e.g., IW itzel et al. 19791; 
IQwen et al. 1983 : ICondon et al. 1984al : ICondonenih*~1989ll . The 

red dashed line marks our threshold to define low-frequency fiat 
spectrum radio sources at ai ow =0A. 

In radio astronomy it is conventional to indi¬ 
cate flat spectrum radio sources as those with a 
two-point spectral ind ex a(u ~0 or typically 
lower than 0.5 (e.g. ICondon et al. 1984al h The ori¬ 
gin of these thresholds reside in the distribution 
of the two-point spectral indices measured between 
~1.4 GHz and ~5 GHz for a number of flux-limited 
so urce sam p les (IWitzel et al. 19791 IQwen et al. 19831 

ICondon et al. 1984;i ICondon et al. 1989H . As shown in 
these analyses, the (unnormalized) spectral-index distri¬ 
butions consist of a narrow steep-spectrum component 
with a{yi,v 2 ) ~0.7 and a broader flat-spectrum com¬ 
ponent centered on a(z/i,i/ 2 ) ~0. As the sample selec¬ 
tion frequency is lowered, the number of steep-spectrum 
sources increases rapidly and the medi an spectr al i ndices 
of bot h co mponents increase (e.g., iKellermann 1974 
ICondon et al. 19891) . As reported by Kellermann et al. 
(1964), the increase in a(iq, 1 ^ 2 ) of each spectral compo¬ 
nent is proportional to the square of its width, so the 
median spectral index of the flat-spectrum component 
changes more rapidly with frequency. 

As shown in Figure 0 even considering three or more 
flux limited subsample of the WSRT-NVSS associated 
sources, we were not able to identify a bimodal behavior 
in the spectral index distribution of ai ow . This, in addi¬ 
tion to the frequency dependence highlighted by Keller- 
rnan (1964), suggest that a different criterion has to be 
chosen to indicate flat spectrum radio sources at low fre¬ 
quencies. 

We noticed that both blazars and Fermi blazars de¬ 
tected in the WENSS show values of ai ow between - 
1 and 0.65 for the largest fraction of their samples 
dMassaro et al. 2013alh Specifically, in our previous 



Fig. 4.— Upper panel) The fractional efficiency g { oLi ow ) defined 
by Eq.[7]for blazars (blue line), the subsample of Fermi blazars (red 
line) all with cti ow between -1 and 0.7 (dashed black line). Lower 
panel) The completeness (j) defined as the ratio between the total 
number of sources with (aiow < a * 0 w ) anc ^ the total number of ex¬ 
pected sources, comp uted for blazars (blue line), the Fermi blazars 
(red line; see Section l4.ll for more details). The chosen threshold 
a* =0.4, adopted to create the LORCAT, is highlighted by the 
dashed vertical line in both panels. 

analysis we considered low-frequency flat spectrum ra¬ 
dio sources as those having ai ow <0.65. This oc¬ 
curred for 90% of the blazars detected by Fermi and for 
mo re than 8 0% of those listed in_the_ROMA-BZCAT 
(jMassaro et al. 20091 iMassaro et al. 201 laf ). However, 
to assemble the LORCAT, we adopted a more conserva¬ 
tive threshold based on the following statistical criterion. 

First we established the number of blazars and Fermi 
blazars that we expect to find within this subsample sim¬ 
ply performing the crossmatch with the ROMA-BZCAT 
within 8".5 as adopted in our previous analysis. We 
found that the number of expected blazars with a WSRT 
counterpart is 979, including 274 Fermi blazars. For a 
given value of the threshold a* ow , we defined the frac¬ 
tional efficiency g(ai ow ) as the ratio between the differ¬ 
ence of total number of sources having: ai ow < a* ow and 
those with (cq 0 w — A a) < a* ow and the total number of 
expected sources N exp within the WSRT-NVSS associa¬ 
tions with -1< ai ow <0.7: 

= N(ai ow < a* ow ) - N{{a low - A a) < a * ow ) 

N exp 

( 6 ) 

where Aq= 0.1. In particular, g(ai ow ) has been com¬ 
puted for all blazars (i.e., gB{oa ow ))i the subsample of 
Fermi blazars (i.e., g 7 (cq OUJ )) with -1< oti ow <0.7 (see 
Figure a. Since the main goal underlying the LOR¬ 
CAT is to have a catalog of potential counterparts for 
the UGSs, we chose the a* ow =0.4 threshold as the value 
corresponding to the peak of the g 7 (cq ow ). According to 
the above threshold the total number of low-frequency 
sources with a flat radio spectrum listed in the LORCAT 
is 28358 having -1< ai ow <0.40. Adopting the above 
criterion on the choice of the cq olu threshold the LORCAT 
catalog will be less complete. However, the selected low 
frequency flat spectrum radio sources are more reliable to 
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be 7 -ray blazar candidates since this criterion ensures to 
avoid the heavy contamination by steep spectrum radio 
sources. 

In Figure [4] we also show the completeness ip of the 
sample considered above defined as the ratio between 
the total number of sources with Our criterion is then 
supported by the comparison at high-frequency described 
in Section [5] ( ai ow < a* ow ) and the total number of 

expected sources: 


Willow) 


N(ui ow < ct low ) 

-V exp 


( 7 ) 


Thus we noticed that for our choice of a* ow =0.4 we 
are able to re-associate 80% of the Fermi blazars of all 
blazars li sted in th e WSRT-NV SS with ai ow between -1 
and 0.65 dMassaro et al. 2013af ). 

In Table Q] we listed all LORCAT sources with their 
WSRT and NVSS names. For all these sources we also 
report the NVSS coordinates, the angular separation R 
between the NVSS and the WSRT positions, the ai ow 
value with its uncertainty ai ow and the WSRT survey 
name which each original source belong to: WISH or 
WENSS. 

Finally, we note that, as f ound by our previous analy¬ 
sis (e.g.' lMassaro et al. 2013ajl . the source density of the 
LORCAT sources is ~1.8 src/deg 2 , given the total 4.7 sr 
of footprint of the combined WENSS-WISH survey (3.1 
sr in the WENSS plus 1.6 sr in the WISH) while accord¬ 
ing to the ROMA-BZCAT the blazar density is currently 
of the order of 0.1 src/deg 2 . So we can expect that only 
about 10% of the sources in the LORCAT are blazar-like, 
however, to confirm this insight optical spectroscopic ob¬ 
servations and high frequency radio data are necessary. 
Moreover, since the ROMA-BZCAT is not a survey and 
it is not a complete catalog, the above estimate on the 
expected fraction of blazars present in the LORCAT has 
to be considered carefully. 


4.2. Flux density distributions at low-frequencies 

Comparing the radio flux densities at 1.4 GHz S '1400 
and at low frequencies Si ow (i.e., 325 for the WENSS 
and 352 MHz for the WISH), as shown in Figure [5j there 
is a good match between the two WSRT and NVSS ob¬ 
servations: bright sources below ~1 GHz tend to be 
among the brightest also above 1 GHz. In Figure [5] 
we also report the line corresponding to a radio spec¬ 
trum of ai ou ,=0. Then in Figure EO we also compare 
the low-frequency radio spectral index ai ow with the 
archival WSRT and NVSS flux densities. The logN- 
logS distribution for all of the WSRT-NVSS associations 
per range of low-frequency spectral indices between -1 
and 1.5 as well as that of our LORCAT are reported in 
Figure [3 Then Figure [ 8 ] shows the logN-logS distribu¬ 
tions for the LORCAT sample for both the WSRT and 
the NVSS flux densities. These logN-logS distributions 
are in agreement with th e evo lution of th e ra dio sourc e 
counts (e.g., iCondon et al. 1984bt (Condon et al. 1998 ). 
These logN-logS distributions computed with both the 
Slow and 5*1400 for all the LORCAT sources appear to 
have the same shape. This is expected as the flux den¬ 
sities are mildly correlated, as shown in Figure [5] In 
Figure [ 8 ] we also show the N oc 5 -15 line expected 
in the case of a uniform source distribution at not too 


WSRT - NVSS correspondences 
LORCAT 


c« 



low 


Fig. 5.— The flux density scatterplot. LORCAT sources (ma¬ 
genta circles) are shown in comparison with those associated in the 
whole WSRT-NVSS crossmatch (black circles). The dashed black 
line marks the radio spectral index oq ow =0. 

large redshift (i.e. in a Euclidean universe). It is well 
know that blazars show a broken luminosity function 
due to the relativistic ef fects of their beamed emission 
(e.g. lUrrv fc Shafer 1984fb this could be also reflected in 
the logN-logS distribution in agreement with that of the 
LORCAT. However to prove this effect redshift estimates 
will be necessary for these low-frequency sources with flat 
radio spectra. 

5. COMPARISON WITH THE GREEN BANK 6-CM RADIO 
SOURCE CATALOG 

A detailed identification of the complete LORCAT 
sample is out of the scope of the present analysis 
and in particular a multifrequency analysis of the op¬ 
tical and the IR counterparts of the LORCAT sources 
will be presented in a sep arate, forthcoming paper 
(jMassaro et al. 2014 in nrep.l ). However to understand 
the nature of the selected low-frequency flat spec¬ 
trum radio sources we performed a crossmatch with 
the G reen B ank 6-cm radio source catalog (GB6)EEI 
(Gregory et al. 1996j l to investigate the spectral proper¬ 
ties of the LORCAT sources at ~5 GHz. 

It is worth noting that among all the radio surveys 
at frequency grater than ~1 GHz the GB 6 is the most 
recent one covering the largest portion of the LORCAT 
footprint since it was performed between 0° and +75° 
in declination. The GB 6 radio source cat alog is also 
complete above 50 mJy (iGregory et al. 19961) . Since the 
CRATES catalog have been compiled using the GB 6 in 
the above range of declination, a comparison with it is 
nested within the following analysis. 

The total number of LORCAT sources within the GB 6 
footprint is 15814. Assuming the difference Aa (i.e., 
Aa = cthigh — oiiow) between the low (i.e., ai ow ) and the 
high-frequency (i.e., cthigh ) spectral indices equal to zero 
with cthigh defined as -1.85dog(5485o/5i4oo), we com¬ 
puted the extrapolated flux density at 4.85 GHz 5 e:ri 485o 
for the LORCAT sources to determine those expected to 
be detected in the GB 6 . We found that above the com¬ 
pleteness threshold of the GB 6 , there are 3219 LORCAT 

16 http://heasarc.gsfc.nasa.gov/W3Browse/all/gb6.html 
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TABLE 1 

LORCAT Main Table (first 10 lines). 


WSRT 

name 

NVSS 

name 

R.A. (NVSS) 

( J2000) 

Dec. (NVSS) 

( J2000) 

R 

arcsec 

a low 

a low 

survey 

0000.0+3323 

J000238+334008 

00:02:38.41 

+33:40:08.3 

6.27 

0.18 

0.12 

wenss 

0000.0+4449 

J000237+450554 

00:02:37.65 

+45:05:54.1 

10.98 

0.38 

0.1 

wenss 

0000.0+5005 

J000236+502220 

00:02:36.82 

+50:22:20.3 

6.78 

-0.28 

0.12 

wenss 

0000.0+6737 

J000235+675422 

00:02:35.79 

+67:54:22.7 

2.78 

0.33 

0.07 

wenss 

0000.0-1838 

J000239-182128 

00:02:39.71 

-18:21:28.8 

29.83 

0.28 

0.16 

wish 

0000.1+4452 

J000244+450928 

00:02:44.17 

+45:09:28.5 

4.39 

0.07 

0.08 

wenss 

0000.1+4628 

J000242+464509 

00:02:42.71 

+46:45:09.0 

10.26 

0.0 

0.09 

wenss 

0000.2-2131 

J000249-211419 

00:02:49.81 

-21:14:19.3 

4.43 

-0.15 

0.03 

wish 

0000.2-2251 

J000250-223437 

00:02:50.77 

-22:34:37.7 

9.52 

0.29 

0.14 

wish 

0000.3+2926 

J000252+294253 

00:02:52.36 

+29:42:53.2 

3.15 

0.22 

0.06 

wenss 


Col. (1) WSRT name. 

Col. (2) NVSS counterpart of the WSRT source. 

Col. (3) R.A. from the NVSS catalog. 

Col. (4) Dec. from the NVSS catalog. 

Col. (5) Angular separation between the WSRT and the NVSS position: R. 

Col. (6) Low frequency radio spectral index oii ow - 

Col. (7) Uncertainty on the o^i ow - 

Col. (8) WSRT original survey: WENSS or WISH. 




Fig. 6. — Left panel) The scatterplot of the low-frequency spectral index ai ow with respect to the WSRT flux density. Right panel) The 
same scatter plot where the NVSS flux density is reported on the x axis. 


sources having Sea, 4850 >50mJy. 

Then searching the correspondences between the LOR¬ 
CAT and the GB6 catalogs, we found 1942 out of the 
3219 (i.e., ~60%) are detected at 6-cm within their posi¬ 
tional uncertainty regions at ler level of confidence, com¬ 
puted between the NVSS and the GB6 positions. In 
particular only 875 out of these 1942 radio sources show 
a flux density S+so above the completeness limit of the 
GB6 survey. The distribution of the high-frequency spec¬ 
tral index ahigh computed with the observed S 48.50 in the 
GB6 for the 2834 LORCAT-GB6 associations is reported 
in Figure|H]together with their Aa histogram. More than 
~75% of the 1942 LORCAT sources detected in the GB6 
still have a “flat” radio spectrum at frequencies above ~1 
GHz (see Figure ©, according to the canonical, widely 
accepted definition of flat spectr um radi o so u rces (i.e., 
oihigh <0.5) (e.g.. IKellermann 1974 iCondon et al. 1989L 
and references therein). This strongly support our def¬ 
inition of low-frequency “flat” radio spectra (see Sec¬ 
tion [XT]). However, a significant fraction (i.e., ~60%)of 
these GB6-LORCAT sources appear to have radio spec¬ 
tra that steepens toward higher frequencies (i.e., Aa > 
°)- 

In particular, it is worth mentioning that the subclass 
of BZQs generally show flat high-frequency spectra, thus 


LORCAT sources with steep high-frequ ency spec tra may 
not actually be BZQs (e.g.. ICondon et al. 198.81 ). How¬ 
ever, Aa > 0 is occurring for a small fraction (i.e., ~5%) 
of the known blazars listed in the ROMA-BZCATand 
these are all classified as BL Lac objects. In Figure [TO] 
we also report the radio spectral index all^ 0 evalu¬ 
ated for all the ROMA-BZCAT blazars that have ra¬ 
dio observations in the NVSS and in the Sydney Univer¬ 
sity Molonglo Sky Survey (SUMSS; iMauch et al. 20031 ) 
at 843 MHz. It is evident that blazars show a clear 
steepening a higher frequencies in agreement with that 
found for the LORCAT sources. In addition, there 
could also be the possibility that t hese radio spectra 
are intrinsical ly mildly cur ved (e.g..[Howard et al 19651 
IKellermann et al. 19691 : iPaulinv-Toth et al. 19721b It is 
known that spectral curvature appe ars at higher frcquen- 
cies i n the_sub millimeter data (e.g.. lGiommi et al. 20071 : 
iGiommi et al. 20121 ). 

Finally, we note the presence of radio sources with 
Aa < 0 in the LORCAT sources might indicate that the 
low-frequency emission could be contaminated by that 
of extended components. These cannot be resolved with 
the large beam of the low-frequen cy survey and in gen¬ 
eral present steep spectra (see e.g JMassaro et al. 201311 
for a recent discussion). 
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Fig. 7. — Left panel ) The logN-logS distributions evaluated for all the WSRT-NVSS correspondences with the low-frequency flux density 
Slow f° r different range of oli ow between the values -1 and 1. Right panel) Same as the left panel for the logN-logS distributions computed 
using the NVSS flux densities at 1.4 GHz. 

LORCAT: logN - log S distributions 




Fig. 8. — The logN-logS distributions of the LORCAT sample 
calculated with the WSRT flux density Si ow (black circles) and 
with that of the NVSS at 1.4 GHz (red squares). A similar shape for 
the two logN-logS distributions of the LORCAT sample is expected 
since the Si ow and 5 i 4 oo flux densities are correlated (see Figurel5|). 
The magenta line indicates the N oc S ~ 1,5 relation expected from 
a uniform source distribution, while the vertical dashed black line 
marks the completeness limit of the WSRT survey at 30 mJy. 

6 . SUMMARY AND CONCLUSIONS 

We have assembled a low-frequency radio catalog of 
flat spectrum sources (LORCAT) built by combining the 
radio observations of the two main WSRT surveys (i.e., 
WENSS and WISH) at 325 MHz and 352 MHz, respec¬ 
tively, with those of the NVSS at 1.4 GHz. The main 
goals underlying the creation of this cat alog are similar 
to those of the CRATES (jHealev et al. 200% since both 
can be used in the future to search for new, unknown 
blazar-like counterparts of the 7 -ray sources EH- 

We defined a new criterion to associate WSRT 
a nd NVSS so ur ces im prov ing our previous analyses 
dMassaro et al. 2013at INori et al. 2014H and we pro¬ 
vided a new definition of flat spectrum radio sources 

17 The LORCAT catalog has been already used for the 7 -ray 
source associations that will be released with the the next Fermi 
catalog actually in preparation. 


Fig. 9. — Left panel) The distribution of the high-frequency spec¬ 
tral index cthigh computed for all the LORCAT sources with a 
counterpart in the GB 6 survey within their radio positional uncer¬ 
tainties (see Section [5] for more details). It is worth noting that a 
significant fraction of sources having flat spectra at low frequen¬ 
cies between ~300 MHz and ~1 GHz (i.e., — 1 < oq ow <0.4) ap¬ 
pear to be relatively flat, according to the general definition (i.e., 
a high <0.5 marked by the vertical black line), also at high frequen¬ 
cies between ~1 GHz and ~5 GHz. Right panel) The distribution 
of the Aq = cx-high ~ a low f° r the LORCAT-GB 6 radio correspon¬ 
dences. Radio sources with flatter high-frequency spectrum with 
respect to the low-energy one have A a < 0 while those steepening 
at high frequencies show Aa > 0. The vertical dashed line marks 
the threshold Aa=0. 

at low-frequencies based on the distribution of the spec¬ 
tral index ai ow between 325 MHz and 1.4 GHz found for 
blazars in the ROMA-BZCAT. Sources with radio anal¬ 
ysis flags as well as double matches between the radio 
surveys have been excluded from our final list. Thus the 
LORCAT sample comprises 28358 radio sources includ¬ 
ing ^667 known blazars having — 1 < ai ow <0.4. 

Then we also compared our LORCAT catalog with the 
the Green Bank 6 -cm (GB 6 ) radio catalog, since it is 
the most recent radio survey covering the largest frac¬ 
tion of the LORCAT footprint at higher frequency (i.e., 
~5 GHz). We found that a significant fraction of the 
LORCAT sources with extrapolated flux densities at ~5 
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BZCAT sources 




Fig. 10.— The distribution of the radio spectral index rvl// 1 for 
all the known blazars that lie within the NVSS and in the SUMSS 
footprints (left panel). The cumulative distribution is shown on 
the right panel where the red dashed line marks the a/// 1 =0.5, 
according to the canonical definition of flat spectrum radio sources. 


GHz above the completeness threshold of the GB 6 are 
detected (i.e., ~ 86 %). In addition they also appear to 
be “flat” spectrum radio sources above ~1 GHz, accord¬ 
ing to t h e cano nical definition (i.e., ahigh <0.5) (e.g., 
iCondon et al. 19891 and references therein). The lack of 
detections for a small fraction of the LORCAT sources 
in the GB 6 footprint could be explained in terms of a 
spectral steeping toward high frequencies (i.e., a mild 
curvature) as already observed in blazars. 

Finally, we highlight that to investigate the nature of 
the LORCAT sources, aiming to identify the fraction of 
7 -ray blazar candidates associable to Fermi sources, a 
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